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Abstract. Smooth linear statistics of random permutation matrices, 
sampled under a general Ewens distribution, exhibit an interesting non- 
universality phenomenon. Though they have bounded variance, their 
fluctuations are asymptotically non-Gaussian but infinitely divisible. 
The fluctuations are asymptotically Gaussian for less smooth linear sta- 
tistics for which the variance diverges. The degree of smoothness is 
measured in terms of the quality of the trapezoidal approximations of 
the integral of the observable. 

1. Introduction 

We study the fluctuations of the spectrum of random permutation matrices, 
or more precisely, of their linear statistics under the Ewens distribution for a 
wide class of functions. The study of linear statistics of the spectrum of ran- 
dom matrices is an active field (for results concerning invariant ensembles, 
see for instance [H], [l6], [l7j, [ISJ, [23], [25j, ^ or ^ and for non- 
invariant ensembles see for instance [6], [7], [9], [11], [19], [26], [27j or [36]). 
All previous results (except [30]) have two common features. Firstly, the 
variance of linear statistics does not diverge for smooth enough functions, 
and thus, no normalization is needed to get a limit law, whereas for less 
smooth functions the variance blows up very slowly (i.e. logarithmically). 
The second feature is that those fluctuations are asymptotically Gaussian 
(except in ^30j again, where invariant ensembles with more than one cut are 
shown to have non-Gaussian fluctuations). We will see that the behavior of 
the variance of the linear statistics of random permutation matrices follow 
the general pattern. But we will also prove that the asymptotic limit law is 
more surprising, in that it is not Gaussian but infinitely divisible when the 
function is smooth enough. This is in contrast to the case where the function 
is less regular, the fluctuations being then indeed asymptotically Gaussian. 
This Gaussian behavior was previously proved by K. Wieand (see }36j) for 
the special case where the linear statistic is the number of eigenvalues in a 
given arc, and for uniformly distributed permutations. 

We first introduce our notations. If N is an integer, Sn will denote the 
symmetric group. We denote by the permutation matrix defined by the 
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permutation a G Sn- For any 1 < i,j < N, the entry Ma-{i,j) is given by 

M^{i,j) = li=^(j). (1.1) 
Mo- is unitary, its eigenvalues belong to the unit circle T. We denote them 

by 

Ai(ct) = e^^'^'^i^'^), . . . , \Nia) = e^*^^^^'") G T, (1.2) 
where ipi{(j), ■ ■ ■ ,ipN{o') are in [0, 1]. 

For any real- valued periodic function / of period 1, we define the linear 
statistic 

N 

I^^nU) := Trf{M^) = ^ f{ipi{a)), (1.3) 

1=1 

where /(e^*'^'^) = f{(p) is a function on the unit circle T. 

We consider random permutation matrices by sampling a under the Ewens 

distribution 

= 9{e + i)...{e + N-iy ^^-^^ 

where 6 > and K{a) is the total number of cycles of the permutation a. 
The case = 1 corresponds to the uniform measure on Sn. 

We study here the asymptotic behavior of the linear statistic Ia,N{f) under 
the Ewens distribution z^tv,^ for any 9 > 0, and a wide class of functions 
/. As mentioned above, the asymptotic behavior depends strongly on the 
smoothness of /. In order to quantify this dependence, we introduce the 
sequence 

RAf)=-Y.f(-) - r/(x)dx. (1.5) 

Using the periodicity of /, it is clear that 

R3U) = -■ (^^/(O) + E / (7) + ^/(l)) - /(^)^^- (1-6) 

So that Rjif) is easily seen to be the error in the composite trapezoidal 
approximation to the integral of / {14] . 

We will see that the asymptotic behavior of the linear statistic /o-,7v(/) is 
controlled by the asymptotic behavior of the i?j(/)'s, when j tends to in- 
finity, i.e. by the quality of the composite trapezoidal approximation to the 
integral of /. The role played by the quality of the trapezoidal approxima- 
tion of / might seem surprising, but it is in fact very natural. It is a simple 
consequence of the fact that the spectrum of the permutation matrix is 
easily expressed in terms of the cycle counts of the random permutation a, 
i.e. the numbers aj{a) of cycles of length j, for 1 < j < A^. Indeed, the 
spectrum of Mq- consists in the union, for 1 < j < A, of the sets of j-th 
roots of unity, each taken with multiplicity aj{a). This gives 

^ ~ ^ ^""^ /k\ 

laMf) = E«^(^) E = E"^(^)E/ - • (1-7) 

j=l u)^=l j=l k=0 ^-^^ 
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So that, using the definition ()1.5p of the Rj^s, and the obvious fact that 
"^j^i jaj{a) = N, it becomes clear that: 

1 N 

Ia,N{f) = N f{x)dx + y^aj{a)jR,{f). (1.8) 

At this point, and using the basic equahty (jl.Sp . it is easy to explain in- 
tuitively the non-universality phenomenon we have uncovered in this work. 
When the function / is smooth enough, the sequence Rj{f) converges fast 
enough to zero to ensure that the linear statistic is well approximated by the 
first terms in the sum (jl.Sp . These terms correspond to the well separated 
eigenvalues associated with small cycles. The discrete effects related to these 
small cycles in the spectrum are then dominant, and are responsible for the 
non-Gaussian behavior. Thus, the appearance of non-universal fiuctuations 
is due to a very drastic localization phenomenon. Indeed, the important 
eigenvalues for the behavior of smooth linear statistics are atypical in the 
sense that they correspond to very localized eigenvectors, those localized on 
small cycles. When the function is less smooth, the variance will diverge 
(slowly) so that a normalization will be necessary. After this normalization, 
the discrete effects will be washed away and the limit law will be Gaussian. 

We will first describe the fluctuations of linear statistics of smooth enough 
functions /, i.e. in the case when the i?j(/)'s decay to fast enough to 
ensure that the variance of the linear statistic stays bounded. 

Theorem 1.1. Let 6 be any positive number, and f be a function of bounded 
variation. Assume that 

oo 

J^ji?,(/)2E (0,oo). (1.9) 
i=i 

Then, 

(1) under the Ewens distribution VNfi, the distribution of the centered 
linear statistic 

i^,N{f)-ni.,N{f)] 

converges weakly, as N goes to infinity, to a non- Gaussian infinitely 
divisible distribution ^f^. 

(2) The distribution fif^g is defined by its Fourier transform 



fijQ{t) = exp j (e^*^ - 1 - itx)dMf{x)^ , (1.10) 
where the Levy measure Mf is given by 

oo ^ 
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(3) The asymptotic behavior of the expectation of the linear statistic is 
given by 

1 N 

nia,N{f)] = N f{x)dx + ey^R,{f) + o{l). (1.12) 

Jo ,=1 

Here, the second term 'Yl!j=iRj{f) 'm-CLy diverge, but not faster than 
logarithmically. 

N 

^i?,(/) = 0(7biiV) (1.13) 
i=i 

(4) The asymptotic behavior of the variance of the linear statistic is given 
by 

N 

Var[I.,;v(/)] = ej;iii,(/)' + o(l). (1.14) 

i=i 

Remark 1. In this theorem (and in the next), we restrict ourselves to the 
class of functions / of bounded variation. This is not at all a necessary 
hypothesis, but it simplifies greatly the statements of the theorems. Our 
proofs give more. We will come back later (in Section 2) to the best possible 
assumptions really needed for each statement. These assumptions involve 
the notion of Cesaro means of fractional order, which we wanted to avoid in 
this introduction. 

Remark 2. We note that the assumption ()1.9p is not satisfied in the trivial 
case where / is in the kernel of the composite trapezoidal rule, i.e when the 
composite trapezoidal rule gives the exact approximation to the integral of 
/ for all j's. In this case, the sequence Rj{f) is identically zero and the 
linear statistic is non-random. Obviously, this is the case for every constant 
function / and for every odd function /, i.e. if 

f{x) = -f{l-x). (1.15) 

It is indeed easy to see then that Rjif) = for all j > 1. 

Remark 3. Consider now the even part of /, i.e. 

/even(:c) = i(/(x) + /(l-x)). (1.16) 

It is clear then that 

RjU) = RjUe.en), (1-17) 

SO that the assumption (|1.9p in fact only deals with the even part of /. 

Remark 4. In order to avoid the possibility mentioned above for all iij(/)'s 
to be zero, we introduce the following assumption 

/even IS not a Constant. (1-18) 

Note that, in general, it is not true that (jl.lSp implies that the sequence of 
i?j(/)'s is not identically zero, even when / is continuous! (See [22] or [28 .) 
But when / is in the Wiener algebra, i.e. when its Fourier series converges 
absolutely, then (jl.lSp does imply that one of the Rj{f)^s is non zero (see 
[28], p. 260). 
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Remark 5. It is in fact easy to compute explicitly the value of the expectation 
and variance of the linear statistic Ia,N{f) for any value of 6 and of N. This 
is done below, in Section [71 The asymptotic analysis is not immediate for 
the values of < 1. 

We now want to show how the assumption ()1.9p can easily be translated 
purely in terms of manageable regularity assumptions on the function / 
itself. 

Corollary 1.2. If f £ , letuj{f',6) be the modulus of continuity of its 
derivative f. Assume that 

oo ^ 

j;^^(/,l/jf <oo, (1.19) 

also assume (|1.18|) in order to avoid the trivial case mentioned above, then 
the conclusions of Theorem hold. 

Of course, the condition (|1.19p is satisfied if / G C^"*"", for < a < 1, i.e. if 
/' is a-Holder continuous. 

We can give variants of the assumptions of smoothness of / given in Corol- 
lary [T21 For instance, 

Corollary 1.3. /// has a derivative in LP, let uj^P\f',6) be the modulus of 
continuity in of its derivative f, i.e. 

u;^PHf',5)= sup j [\f'{x + h)-f'ix)A '\ (1.20) 

Assume that 

Jp\f',6)<6" witha>-, (1.21) 

P 

also assume (jl.lSp in order to avoid the trivial case mentioned above, then 
the conclusions of Theorem hold. 

It is of course also possible to relate the i?j(/)'s to the Fourier coefficients 
of /. Indeed, if the Fourier series of / 

oo oo 

f{x) = oq + a„ cos(n27rix) + 6„ sin(n27ri2;) (1.22) 

n=l n=l 

converges, then the Poisson summation formula shows that 

oo 
n=l 

Using this relation, it is easy to prove the following Corollary: 

Corollary 1.4. // / is in the Sobolev space , for s > 1, and if one 
assumes (jl.lSp . then the conclusions of Theorem \l.l\ hold. 

Remark 6. The formula 11.231 gives an expression for the asymptotic variance 
of the linear statistic 

oo oo 

hm Var[I,,jv(/)] = ^ Vji?,(/)' = 9 V akaid{k,l), (1.24) 

N—^oo ^ — ' ^ — ' 
j=l k,l=l 

where d{k, I) is the sum of the divisors of the integers k and I. 
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We now give two interesting examples of functions satisfying the conditions 
of Theorem 11.11 



Example 1. Let / be a trigonometric polynomial of degree k. Then, Rj{f) = 
for all j > k. Obviously, the condition (jl.9p of Theorem 11.11 is satisfied 
and the limit distribution fij^Q is a compound Poisson distribution with Mj 
given by 

k , 



Example 2. Let / G C°° and / = 1 on [a, b] and / = on [a — e,b + e]^, then 
the result of Theorem 1 1 . 1 1 applies . So, the centered linear statistic Iaj^{f) — 
E[/o-^(/)] has a finite variance and a non-Gaussian infinitely divisible limit 
distribution. This is a very different behavior from the case / = l[a,fe] (see 
below), where the limit is Gaussian. 

We now give our second main result, i.e. sufficient conditions ensuring that 
the variance of the linear statistic Ia,Nif) diverges and that the linear statis- 
tic converges in distribution to a Gaussian, when centered and normalized. 

Theorem 1.5. Let 9 be any positive number, and f be a bounded variation 
function such that 

oo 

Y^jRjif? = oo. (1.26) 

i=i 

Then, 

(1) under the Ewens distribution VNfij the distribution of the centered 
and normalized linear statistic 

J.,^(/)-E[/^,^(/)] 

^Marl^^Nif) 

converges weakly, as N goes to infinity, to the Gaussian standard 
distribution AA(0, 1). 

(2) The asymptotic behavior of the expectation of the linear statistic is 
given by 

1 N 

ni^,N{f)] = N f[x)dx + ey^R.j{f) + 0{l). (1.28) 

Here, the second term ^^i-Rj(/) may diverge, but not faster than 
logarithmically. 

N 

Y,RAf) = 0{\ogN) (1.29) 
j=i 

(3) The asymptotic behavior of the variance of the linear statistic is given 
by 

N 

Var[/,,^(/)] ^eJ2jRjiff- (1-30) 



EIGENVALUES OF RANDOM PERMUTATION MATRICES 



7 



Example 3. Consider / = l(a,b) for an interval (a, 6) C [0,1]. Ia,N{f) is 
then simply the number of eigenvalues in the arc [e^*™, e^*'^^]. The function 
/ is obviously of bounded variation. This example has been treated in the 
simple case where 9 = 1 in [3B] . We will see here that Theorem 11.51 enables 
us to extend the results of [36] to any value of > 0. Indeed, the error in 
the composite trapezoidal approximation Rj{f) is very easy to compute for 
an indicator function: 

Rj{f) = ]{{ja}-{jb}). (1.31) 
Obviously, in this case. Theorem 11.51 applies and we have that 

VC9 log N ^-'nu,!;. 

We can also deduce the asymptotic behavior of the expectation and of the 
variance, using the conclusions of Theorem 11.51 and the computations made 
in the particular case ^ = 1 in [36]. Indeed, it is shown in [36j, that for a 
constant ci(a, b) 

N 

Y^Rjif) = -ci{a,b)logN + o{\ogN). (1.32) 

So that from the statement proven in Theorem 11.51 : 

1 N 

ni^,Nif)] = N f{x)dx + 9y^R,{f) + 0{l). (1.33) 

We see that 

WaAf)] = -a)- 9ci log N + o(log N). (1.34) 

The value of ci(a, b) is studied in [36]. It depends on the fact that a and b 
are rational or not. It vanishes if a and b are both irrational. 
We also have, from the computations in [36], that there exists a positive 
constant C2(a, 6) such that 

N 

Y^jRjUf = C2{a, b) log N + o(log N). (1.35) 
i=i 

So that we have for any ^ > 0, by Theorem II. 5| that 

Var/,,,v(/)~C2(a,6)01ogiV. 

The value of C2(a, 6) also depends on the arithmetic properties of a and b, 
and is studied in [36]. 



Remark 7. We want to point out that / being of bounded variation is not a 
necessary condition in order to get a Gaussian limit distribution. But, when 
/ is of bounded variation, it is easy to see that there exists a constant C 
such that 

Varl,,jv(/) < ClogiV. (1.36) 

The case treated in the example above gives the maximal normalization for 
functions of bounded variation. 
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The remainder of this article is organized as fohows. In Section 2, we state 
our results with weaker assumptions than the theorems given in this in- 
troduction. These assumptions use the classical notion of Cesaro means 
of fractional order, which we recall in the first subsection of Section 2. In 
Section 3, we prove the Corollaries 11.21 11.31 and 11.41 using estimates on the 
trapezoidal approximation. In order to prove the main results of Section 2, 
our main tool will be the Feller coupling. This is natural since the problem 
is translated by the basic equality (jl.Sp in terms of cycle counts of random 
permutations. In Section 4, we will need to improve on the known bounds 
for the approximation given by this coupling (see for example [3] or [8]) 
and relate these bounds to Cesaro means. We will then be ready to prove 
in Section 5 and Section 6 our general results as stated in Section 2 and 
that these more general results imply the two theorems of this introduction 
Theorems 11.11 and 11.51 . Finally in the very short Section 7, we give an ex- 
plicit expression for the expectation and variance of the linear statistics as 
promised in Remark [5j 

2. Cesaro means and convergence of linear statistics 
2.1. Cesaro Means. 

We will state here our optimal results in terms of convergence of the Cesaro 
means of fractional order. First, we will need to recall the classical notion 
of Cesaro means of order 9 and of Cesaro convergence (C, 6) for a sequence 
of real numbers, say s = (sj)j>o (see [38], Volume 1, p. 77, formulae (1-14) 
and (1.15)). 

Definition 2.1. (i) The Cesaro numbers of order a > —1 are given by 

(ii) The Cesaro mean of order 6* > of the sequence s = {sj)j>o is given 

by 

<'Nis) = E^'r (2.2) 

(iii) A sequence of real numbers s = (sj)j>o is said to be convergent in 
Cesaro sense of order 9 (or in (C, 6) sense) to a limit £ iff the sequence 
of Cesaro means o"^(s) converges to i. 

Let us recall the following basic facts about Cesaro convergence (see [55]): 

Lemma 2.2. (i) Convergence in the {C,6i) sense to a limit i, implies 
convergence (C, ^2) to the same limit for any 61 < 62- 

(ii) Usual convergence is (C,0) convergence. The classical Cesaro con- 
vergence is (C,l) convergence. 

(iii) If the sequence {sj)j>o is bounded and converges {C,Oi) to a limit i 
for some value 9i > 0, then it converges (C, 9) to the same limit, for 
any 9 > 0. 

These facts are all classical, see [38] for a proof, in particular Lemma (2.27), 
p. 70, Volume 2 for a proof of (iii). 
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2.2. The case of bounded variance, non-Gaussian limits. 

We will give here a sharper statement than Theorem 11.11 and prove that it 
implies Theorem ll.il Define the sequence u{f) = {uj{f))j>i = {jRj{f))j>i. 

Theorem 2.3. Let 9 be any positive number, and assume that the sequence 
l^(/)l = (l%(/)l)i>i converges to zero in the Cesaro {C,9) sense if 6 < 1. 
Also assume that 

oo 

Y,jRjif? ^iO,oo). (2.3) 
i=i 

Then, 

(1) under the Ewens distribution i^N,e, the distribution of the centered 
linear statistic 

I^,N{f)-nia,N{f)] 

converges weakly, as N goes to infinity, to a non- Gaussian infinitely 
divisible distribution fiffi- 

(2) The distribution fif^g is defined by its Fourier transform 

fij Q{t) = exp j (e^*^ - 1 - itx)dMf{x)^ , (2.4) 
where the Levy measure Mf is given by 

oo ^ 

i=i 

(3) The asymptotic behavior of the expectation of the linear statistic is 
given by 

1 N 

WMf)] = N / f{x)dx + V Rjif) + o(l). (2.6) 

Here, the second term X^jLi if ) fno-U diverge, but not faster than 
logarithmically. 

N 

Y,Rjif) = 0{V^) (2.7) 
j=i 

(4) If, on top of the preceding assumptions, one assumes that the se- 
quence u{f)'^ = {uj{fY)j>i converges in Cesaro (C, 1 /\ 6) sense, 
then the asymptotic behavior of the variance of the linear statistic is 
given by 

N 

Var[/.,^(/)] = 0j]ji?,(/)2 + o(l). (2.8) 
Theorem 12.31 will be proved in Section 5. 
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2.3. The case of unbounded variance, Gaussian limits. 

We will give here a slightly sharper statement than Theorem 11.51 and prove 
that it implies Theorem 11.51 

Theorem 2.4. Let he any positive number, and assume that 

oo 

^ji?,(/)2 = oo (2.9) 



and that 



^ma^lj-Rjl =o(??iv), (2.10) 



where rj% = eZf=ijRj if?- Then, 

(1) under the Ewens distribution VNfij ihe distribution of the centered 
and normalized linear statistic 

h,N{f)-Wa,N{f)] 



VVar/,,Ar(/) 



2.11 



converges weakly, as N goes to infinity, to the Gaussian standard 
distribution A/'(0, 1). 
(2) The asymptotic behavior of the expectation of the linear statistic is 
given by 



W<^,N{f)] = N / f{x)dx + y^R,{f) + o{'nN). 

Jo ,=1 



N 

2.12) 



Here, the second term X^jLi-^i(/) diverge, but not faster than 
logarithmically. 

N 

Y,Rj{f) = o{vNVl^) (2.13) 
i=i 

(3) The asymptotic behavior of the variance of the linear statistic is given 
by 

N 

Y^T[h,N{f)] ~ v% = eY,jRjiff- (2-14) 

This theorem will be proved in Section 6. 



3. Estimates on the trapezoidal rule and proofs of the 
Corollaries II. 2L 11.31 and 11.41 

In this section, we will discuss known results about the quality of the com- 
posite trapezoidal approximation for periodic functions, in order to relate 
the decay of the i?j(/)'s to the regularity of /. Moreover, we will give proofs 
of Corollary [Ol Corollary O Corollary (Ll 
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3.1. Jackson-type estimates on the composite trapezoidal approx- 
imation. 

In order to relate the decay of the i?j(/)'s to the regularity of /, we can 
use two related approaches. First, we can control directly the size of the 
i?j's by Jackson type inequalities as in [10], [13] or [31]. Or we may use the 
Poisson summation formula given in ()1.23p and use the decay of the Fourier 
coefficients of /. 

We start by using the first approach, and recall known Jackson-type es- 
timates of the error in the trapezoidal approximation. 

Lemma 3.1. (i) There exists a constant C < 179/180 such that 

\Rj{f)\<Cu2{fM2j), (3.1) 

where 

L02{f,5)= sup \f{x + 2h)-2f{x + h) + f{x)\. (3.2) 

\h\<S,x£[0,l] 

(ii) // the function f is in , then 

\Rm < C ""^^' ■ (3-3) 

(iii) // the function f is in W^'P, then 

|i?,(/)|<0)(/M/i)-^. (3.4) 

Proof. The first item is well known, see (see jlOj). 

The second item is a consequence of the first, since by the Mean Value 
Theorem 

u;2{f,5) <6oj{f',26). (3.5) 
The third item is also an easy consequence of the first since 

f{x + 2h)-2f{x + h) + f{x)= {f{t + h)-f'{t))dt. (3.6) 

J X 

So that 

\f{x + 2h) - 2f{x + h) + f{x)\ < (^J^ \f'{t + h)- f'{t)\Pdt^ ' h"-ir, (3.7) 
which shows that 

a;2(/,5) <>)(/', 5)5'~'/P. (3.8) 

□ 
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3.2. Proofs of Corollary 11.21 and 11.31 using Jackson bounds. 

Proof of Corollarv M.^ We can control the decay of the iij(/)'s using the 
item (ii) of Lemma 13. H which implies that 

jR,{ff < C'^^^i^^ (3.9) 

It is then clear that under the assumption (|1.19|) . the series J^JLiJ^jif)'^ 
is convergent. But p.l9p implies that the Fourier series of / is absolutely 
convergent, so by the result mentioned above ([2H], p. 260) it is true that 
p.l8|) implies that one of the i?j(/)'s is non zero. And thus, J^'jLijRjif)'^ ^ 
(0,oo). If we add that / is obviously of bounded variation, we have then 
checked the assumptions of Theorem 1 1 . 1 1 and thus, proved Corollarv ll.2l □ 

Proof of Corollary \1.3[ We can here control the decay of the iij(/)'s using 
the item (iii) of Lemma |3.H and the assumption (jl.2ip . which imply that 

So, if Q > 1/p, the series Yl'jLi jRjif )'^ is convergent, since 

^•^|(/) ^ ,i+2(li/p) • (3-11) 

Moreover, as above, it is easy to see that (|1.2ip implies that the Fourier series 
of / is absolutely convergent, so by the result mentioned above ([28], p. 260) 
it is true that (jl.lSp implies that one of the i?j(/)'s is non zero. Again, / 
is obviously of bounded variation, we have then checked the assumptions of 
Theorem 11.11 and thus, proved Corollary 11.31 □ 

Remark: It is in fact true that linij^oo jPj if) = is satisfied as soon as 
/ G W^^P (see [ISj). 

3.3. Proofs of Corollary 11.41 and the Poisson summation formula. 

We now turn to the proof of Corollary 11.41 using the second possible ap- 
proach, i.e. the Poisson Summation Formula, ()1.23p . 

Proof of Corollary \1.4\ Let / be in H^, s > 1 and consider its Fourier series 

oo oo 

f{x) = ao + On cos(n27rix) + 6„ sin(n27rzx). (3.12) 

n=l n=l 

Then there exists a sequence {ck)k>i S such that 
So, 



Qk = (3.13) 



e>i * 

is in £^ by Lemma 4 of [31], p. 131. Thus, using the Poisson summation 
formula (fOHD . 



Rjif) = , (3.15) 



J' 
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which is more than enough to prove that the series Yl'jLij^jif)'^ is con- 
vergent. Moreover, as above, it is easy to see that (jl.lSp imphes that one 
of the -Rj(/)'s is non zero, and that / is obviously of bounded variation. 
We have then checked the assumptions of Theorem 11.11 and thus, proved 
Corollary Ol □ 

4. Bounds on the Feller coupling and Cesaro Means 
4.1. The Feller Coupling. 

Let cr G Sj\f be a given permutation and aj{a) be the number of j-cycles 
of a. A classical result is that under the Ewens distribution the joint 
distribution of {ai{a), . . . , ajyio')) is given by 

N 



J^7V,0[(ai(<T),---,aAf(o-)) = {ai,...,aN)] = lE'Lii<*.=A^fl~ 11 (t) 77' 

(4.1) 



where Of^N) = 0{9 + 1) . . . {9 + N - 1). 



We recall now the definition and some properties of the Feller coupling, a 
very useful tool to study the asymptotic behavior of aj (a) (see for example 
m, p. 523). 

Consider a probability space {Q,T,F) and a sequence (^j)j>i of independent 
Bernoulli random variables defined on {0,,J^) such that 

m = 1] = . , ^ , and mi = 0] '"^ 



+ ' 9 + i-l 

For 1 < j < A^, denote the number of spacings of length j in the sequence 
16---Civl hy CjiN), i.e. 

N-j 

CjiN) = ^ . . . il-C^+J^l)C^+jHN-j+lil-^N-J+2) • • • (I-^n)- 

1=1 

(4.2) 

Define {Wjm)j>i by 

oo 

Wjm= - Ci+l) ■ ■ ■ {I - (4.3) 

i=m+l 

and set for j > 1, 

W, := Wjo. (4.4) 

Define 

Jn = min{j > 1 : (,N-j+i = 1} (4.5) 

and 

Kn = mm{j > 1 : ^n+j = !}• (4.6) 
With the notations above, we state the following result of [8j, p. 169: 
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Theorem 4.1. Under the Ewens distribution vn^j 

(i) {Cj{N))i<_j<_N has the same distribution as {aj{a))i<j<N , i-^- for 
any a = (oi, . . . , a^) G N^, 

¥[{Ci{N), . . . ,Cn{N) = a] = UN,eKai{^), ■ ■ ■ ,aN{cT) = a], (4.7) 

(ii) {Wj)i<j<N CLf^ independent Poisson random variables with mean 
0/j, 

(iii) and 

\Cj{N) - Wj\ < WjN + l{j,+x^=,+i} + l{j,=,}. (4.8) 

We will need to improve on the known results for the Feller coupling. In 
particular we will need the following. For any sequence of real numbers 
{uj)j>i, define 

N 

GN = Y,^jCj{N) (4.9) 

i=i 

and 

N 

i=i 

We will need to control the and L^-distances between the random vari- 
ables G]\f and Hi\f. In order to prove Theorem 11.11 and Theorem 11.51 we 
will apply these estimates to the case where the sequence Uj is chosen to be 
u,U)=jR,{f). 

4.2. bounds on the Feller Coupling. 

We begin with the control of the L^-distance in this subsection. We first 
state our result in a very simple (but not optimal) shape. 

Lemma 4.2. For every 9 > 0, there exists a constant C{9) such that, for 
every integer N , 

- Hn\) <Cie) laa.^ \uA. (4.11) 

l<j<N ■' 

This result is a trivial consequence of a deeper result, that we now give after 
introducing some needed notations. We recall that for any real number x 
and integer fc, 

'x\^ ^ x{x-l)...{x-k + 1) 
We now define for any > and every 1 < j < A^, 

*^^^)^=V N-j )[ N ) =^^ e + N-k-r ('-^'^ 

where 7 = — 1. 
We then have: 

Lemma 4.3. 

-Hn\<^Y. I^^-I + nY1 l^.-l^^(i) (4-14) 
i=i i=i 



E|C,-(iV) -Wj\<^ + -^n{3) (4.15) 
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Lemma 14.31 is obviously a direct consequence of the following: 

Lemma 4.4. Let 6 > 0, there exists a constant C{9), such that, for every 
l<j<N 

N ' N 

In order to prove Lemma |131 we note that, bv \4:.8\ 

E\Cj{N) -Wj\< E{WjN) + IP(JiV + KN=j + l)+ P(Jv = j). (4.16) 
It thus suffices to provide bounds on E[WjAr], f[J]\f = j] and F[Jj\f + Kjy = 

Lemma 4.5. For any 9 > and for every i < j < N , 

nWjN) < (4.17) 



Proof. Let 
then, for i > 2, 



' := &(1 - £i+i) ^ ^ ^ (1 - Hi+i-iK.+l, (4^18) 



By (j4.19p . we have immediately that, for any 6 > 0, 

i=N+l £=N 

□ 

We compute next the distribution of the random variable Jn explicitly. 
Lemma 4.6. 

njN=j] = ^^N{j). (4.21) 

Proof. The random variable Jn is equal to j if and only if = 0, ^Ar_i = 
0, . . . , ^N-j+2 = and Cv-j+i = 1- So, for any 1 < j < iV, 

N-l N-2 N-(j-l) ft 



+ iV-l 6l + A^-2'"6l + A^-(j-l) e + N -j 

fc=0 

which proves the claim. □ 
We now bound the distribution of the random variable Jm + Kn- 
Lemma 4.7. For any 9 > 0, 

nKN + JN=j + l]<-^. (4.23) 
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Proof. Consider first the random variable K^. For any ^ > 0, 



< 



< 



N N + 1 

+ N ^ 9 + N + 1 "' 6 + N + j - 2 " e + N + j -I 
NO 



' + N + j-2){e + N + j-l) 



< 



+ N + j -I - N + 9 
For any 6* > 0, use (|4.24p to write 

i 

P[KjV + JiV = j + 1] = Y.^{JN=j + l-mKN=^) 



(4.24) 



£=1 



< 



N + e 



Y,njN=j + i-i) 



1=1 



(4.25) 
□ 



The last three lemmas imply the result of Lemma 14. 4[ We have now con- 
trolled the L^-distance between Gat and H^. 

4.3. bounds on the Feller coupling. 

We now turn to the control of the L^-distance between the random variables 
Gm and Hj^j. We first state our result in a simple (but not optimal) shape. 

Lemma 4.8. For every 9 > 0, there exists a constant C{9) such that, for 
every integer N , 



E {{Gn - Hn?) < G{9) max 



(4.26) 



i<i<Af 

This result is an immediate consequence of the following much more precise 
statement. 

Lemma 4.9. For every 9 > 0, there exists a constant G{9) such that, for 
every integer N, 



E{{GN-HMf) < G{9) 



N N 



N 



N 



N 



j=i k=i 
1 ^ 

-^|n,f*^(j) 



Proof. We note that 



N 



¥.{{GN-HNf)< \uj\\uk\n\Cj-W,\\Gk-Wk\). 
j,k=i 



(4.27) 



(4.28) 
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By (gSj), for any fixed l<j,k<N, 

\Cj - Wj\\Ck - Wk\ < Wj,NWk,N + Wj,Nljr,=k + Wj,Nlj^+K,,=k+l 

+ '^JN+KN=j + l^k,N + tjN+KN=j + l^JN=k 

+ 1 jjv+_ffjv=i+i 1 jjv+i<:jv=fc+i 

+lj^=jlj^+K^=k+i. (4.29) 

To control ()4.29p . we will give upper bounds for all the terms on the RHS. 
We start by giving a bound for E{Wj^]\[Wk,N)- By (|4.18p . we have 

i,l>N+l 

i,i>N+l i,e>N+l i>N+l 

(4.30) 

We write the first term on the RHS as follows: 

E f^^^'^^f + E f^^^'^^f + E U^'^U^'?r (4-31) 
i,e>N+l i,e>N+l i>N+l 

i<e<i+j i+j<e 

It is easy to see that for any £ G (i, i + j), U^^'^uf^'* = 0. If £ is strictly larger 
than i+j, then u''/^ and uf^^ are independent. This gives, using ()4.19p . 



i,£>N+l i,l>N+l 
i+j<l i+j<i 



<- E^(^i")E^ 

^ 1 cie) , 



Also, by the same argument, 



E f^S- < E iE(eoiE(e.+,)iE(e^+,+fc 

i=Af+l i=Af+l 

- . 4:^ . -l){i + j + k-l) - ^ ^ 



oo 



=N+1 

For the bound of (|4.3Up . we consider now the second and the third term on 
the RHS. But the second term can be bounded similarly to the first term. 
For the third term in ()4.30p . we observe that U^'^^uf''^ = ii j ^ k. So, by 

(iron 
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>Af+l i>Af+l ^ ^ i>Ar+l 

g ^ c{e) 

(i- 1)2 - AT 



i>Af+l 

This gives 



So, 



AT / 1 ^ \ 1 ^ 

Y,Wj\H\nwj,NWk,N) < ci{e)\-Y.\''A +^2W^E 

i,k=l \ j=l J j=l 

(4.36) 

Obviously, Wj^n and lj^=fc are independent. So, the expectation of the 
second term on the RHS in (j4.29p is bounded as fohows: 

E{Wj,Nlj^=k) < ^nJN = k). (4.37) 
Of course, this bound is also valid for E(l j^=jVFfc^iv). 

Then, 

N N N 

j,k=l j=l k=l 

(4.38) 

For Wj^NtjN+KN=k+i, we write 

E(PFj,Arlj^+Xjv=fc+i) =E ^E ^j,N'^JM+KN=k+i'^JN=e^ 
k 

= E nWj,N+k~elK^=k+i-e)nJN = e). (4.39) 



But, 

=k+i-e) 



E ( E U^'^^K^=k+i-i] + E [ul^l,^,_,lK^=k+i-e) , (4.40) 

\i>Af+fe+l-^ / 



where Yli>N+k+i- 

_,UI'^ and lK^=k+ i_f are independent and 

= E (.^Ar+fc+i_^(l - iN+k+2~d ... (1 - iN+k+j-l)iN+k+j+l-iiN+k+l-i) 

- ^ (^^v+fc+i-^ j ^ (AT + A; - ^)2 - N^- ^^-^^^ 
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So, by ^M) 

E{W,,N+k-elK^=k+i-i) < ^F{K^ = k + l-£) + ^<^, (4.42) 



which gives 



E{Wj^N+k-etK^=k+i^i) < (4.43) 



This gives also the bound for K{lj^+KN=j+i^k,N)- 
Then, 

Wj\WkmWj,Nlj^+Kr,=k+i) < Cie)l-^\u,\ \ .(4.44) 
For the remaining terms in (|4.29|) . we observe that 

E{lj^+K„=j+lljr,=k) = E{lK^=j+i-klj^=k)=nKN=j + l-k)¥{JN = k) 

< ^PiJ^ = k). (4.45) 
This apphes for ljj^=jljj^+KN=k+i, as well. 
Then 

N N N 

j,k=l j=l k=l 

(4.46) 



Also, 



so, 



( 1 a k = j I A 

ljjv+ii-jv=j+lljjv+/fjv=A;+l = I Q if^T^j, ' 



and by (iOajl 

N N 

Y Wj\\uk\E{ljj^+KN=j+l'^JN+KN=k+l) = ^ |'UjpP(JAr + Kat = j + 1) 

j,k=l j=l 

< ^El-/. (4.49) 



N 



It is obvious that lj^=jlj^=fc = for /c 7^ j. So, 

E,l.„.,..„..) ^ I (4.0, 
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Then, 

N N 

Y,WjMJN=j), (4.51) 

j,k=l j=l 

which, using also Lemma 14.61 proves the claim of Lemma 14.91 □ 
4.4. Cesaro Means and the Feller Coupling Bounds. 

The hnk between our estimates and Cesaro means of fractional order is given 
by an interesting interpretation of Cesaro means of order 9 in terms of the 
random variable Jat. 

Lemma 4.10. The Cesaro mean of order 6 of a sequence s = {sj)j>o, 
with So = 0, is given by 

N ^ 6 ^ 

i=i i=i 



The proof of this lemma is immediate from Lemma 14.61 the Definition 12.11 
of the Cesaro means and of the numbers ^ given in ()4.13p . 
Using this interpretation of the Cesaro means, we can state our results 
about the and distance between the variables Gn and Hn given 
in Lemma 14.31 and Lemma 14.91 in terms of the Cesaro means of the sequence 
Uj{f) and Uj{ff. 

Theorem 4.11. For an 9 > 0, there exists a constant C{9) such that 
(i) 

E|G^ - Hn\ < C{9){a],{\u\) + a%{\u\)) (4.53) 

(ii) and 

E {{Gm - Hr,f) < C{9)[<y],{\u\f + a],{u') + ^^(|n|)a^(|^|) + a^iu")]. 

(4.54) 

This theorem is simply a rewriting of Lemma [4. 3 1 and Lemma [4. 9 1 using the 
identification given in Lemma 14.101 It implies easily the following results 

Theorem 4.12. If the sequence {uj)j>i converges in Cesaro {G,9Al) sense 
to 0, then 

lim ^Gn -Hn\ = 0. (4.55) 

TV— >-oo 

Proof. By assumption, the sequence converges in (C, 1) and in (C, 9) sense 
to 0. Thus, the RHS of the bound given in Theorem 14.111 tends to zero, 
which proves Theorem 14.121 

□ 



Similarly we can get the following result about convergence in L 

-tj|)j>i and {u'j)j 



Theorem 4.13. If the sequences (|%|)j>i and (u^)j>i both converge to zero 



in Cesaro (C, 9 then 



lim E((Gjv - HN f) = (4.56) 
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Proof. By assumption, the sequences (|tij|)j>i and {u'j)j>i converge in (C, 1) 
and in (C, 0) sense to 0. Thus, the RHS of the bound given in (|4.54p tends 
to zero, which proves Theorem 14.131 □ 

5. Proof of Theorem 12.31 and Theorem 11.11 

5.1. A simple convergence result for series of Poisson random vari- 
ables. 



We give here a result of convergence in distribution for the random variables 

N 

HN{f) = Y.W,u,{f), (5.1) 
i=i 

to an infinitely divisible law. This result is elementary since it only uses the 
fact that the random variables Wj's are independent and Poisson. 

Lemma 5.1. Under the assumption (|2.3p . i.e 



N 

G (0,oo), (5.2) 

i=i 

the distribution fiN of Hj^j —¥j[Hjsi] converges weakly to the distribution fif g 
defined by ([23]). 

Proof. The Fourier transform of Hjsf — K[Hiy] is easy to compute, indeed: 



log/lAr(t) = logE 



,it(i/jv-E[i/jv]) 



N 



logJjE 



exp ituj Wj — : 



N 



Y^l^e^'-^j - ituj - 1). 
i=i 



Obviously, for \t\ < T, 

(e**«. - itUj - 1) 



< 



j 2 



2 j 



Bv ll.91 log 'flN{t) converges absolutely uniformly and its limit 



3 



(5.3) 



(5.4) 



(5.5) 



is continuous. By Levy's Theorem, exp(?/^(t)) is the Fourier transform of the 
probability measure /if^g and ij,n converges in distribution to /i/^g as N goes 
to infinity. □ 

Obviously, /x j^g is an infinitely divisible distribution and its Levy-Khintchine 
representation is easy to write. We recall that an infinitely divisible distribu- 
tion /X has Levy-Khintchine representation (a, M, cj^) if its Fourier transform 
is given by 

Jl{t) = exp f / f e**^ - 1 - dM{x) + iat - \a'^tA , (5.6) 
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where a G M, cr > and M is an admissible Levy measure, i.e. 

r 

/ KdM(x) < oo. 

J l + x^ 

The distribution fif^g in Lemma l5. II has therefore a Levy-Khintchine repre- 
sentation (a, 9M, 0) with 



and 



oo ^ 



It is easy to see that the assumption 11.91 imphes that J x^dM{x) < oo so 
that M is admissible. 

5.2. Proof of Theorem 12731 

We proceed now to the proof of Theorem 12.31 by using the Feller coupling 
bounds proved in Section 4. 

We first prove the first and second statements of Theorem 12.31 Under the 
assumption that X^j^i j-Rj < oo, we have seen that the sequence \uj\ con- 
verges in (C, 1) sense to zero. Moreover, if < 1, we assumed in Theorem l2.3l 
that the sequence \uj\ converges in {C,6) sense to zero. Thus, we know that 
the assumption of Theorem 14.121 is satisfied, and thus that, 

lim E\Gn-Hn\ = 0. (5.9) 

Using now Lemma l5. 11 we have proved that Gn — E{Gn) converges in dis- 
tribution to ^ffi defined by ()2.4p . But, by the basic identity (jl.Sp . we know 
that Ia,N{f) — ]E[/(j,iv(/)] has the same distribution as Gjv(/) — E[GAr(/)]. 
This proves the first two statements of Theorem 12.31 



The third statement is simple. Indeed, by (jl.Sp and by (|5.9p . 

nia,N{f)]=N f' f{x)dx+E{GN) = N C f{x)dx+E{HN) + o{l). (5.10) 
Jo Jo 

In order to complete the proof, it suffices to mention that the expectation 
of i^AT is easy to compute: 

N N 

i=i i=i 
This proves the third statement of Theorem [ 



The proof of the fourth statement follows a similar pattern. Again, by (jl.Sp . 

Var[/,,,v(/)] =Var(G^). (5.12) 
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But if one also assumes, as in the fourth item of Theorem 12.31 that the 
sequence (uj) converges in (C, 1 A 9) sense to zero, then by 14.131 we know 
that 

lim E{{Gn - Hn?) = 0. (5.13) 

iV->-oo 

This, and 15.91 imply that 

Var(G7v) = YaT{HN) + o(l) (5.14) 
In order to complete the proof, it suffices to compute the variance of Hj\[: 



Yax{HN) = e^^ = e^jR] (5.15) 
This proves the fourth statement and completes the proof of of Theorem 



5.3. Proof of Theorem 11.11 

We show here how Theorem 12.31 implies Theorem II. li 
We will need the following simple facts. 

Lemma 5.2. (i) The assumption X^j^i i-Rj(/)^ < oo implies the (C, 1) 
convergence of the sequence {\uj{f )\)j->i to zero. 

(ii) // one assumes that Yl'jLij^jif)^ < and that the function f is of 
bounded variation, then the sequence {\uj{f)\)j>i converges in {C,6) 
sense to zero, for any 9 > 0. 

(iii) // one assumes that X^^i J-Rj(/)^ < oo and that the function f is of 
bounded variation, then the sequence {uj{f)'^)j>i converges in {C,6) 
sense to zero, for any 9 > 0. 

Proof. The first item is well known (see statement (a), p. 79 of [38j, Volume 
1). It is a consequence of the simple application of the Cauchy-Schwarz 
inequality 

TV 1 ^ ^ °° 



'N ^ " - N 

j=i j=i j=i j=i 



So that 



limsupl-^^/,1 < (5^ii?|)i (5.17) 

But the LHS of (|5.17|) does not depend on the initial k values of the sequence 
Uj. By setting these k values to zero, and by taking k large enough, we can 
then make the RHS as small as we want. This implies that 

1 ^ 

lim — y = 0. (5.18) 

This is the (C, 1) convergence to zero, claimed in item (i). 

In order to prove the item(ii), we need the following observation. 
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Lemma 5.3. // the function f is of bounded variation, then 

\Rm < (5.19) 

where TV{f) denotes the total variation of f . 

Proof. Since / is of bounded variation, it can be written as a difference of 
two non-decreasing functions 

/ = /+-/-. (5.20) 



^^•(/^) = E / ' (/^(-) - f^^^))d^- (5-21) 
7 — f^J- 3 



Using ()1.5 



So that 



k=0 



3 



\RAn\ < - E - /^(-)) < V(/+). (5.22) 

Using the same argument for gives the result of Lemma l5.3i □ 

So, this shows that the sequence {uj{f))j>i is bounded, when / is of bounded 
variation. Now, using the item (i) of this Lemma and item (iii) of Lemma[221 
we see that the sequence (|%(/)|)j>i and thus, Uj{f) converges in {C,6) 
sense to zero, for any value of > 0. 

The last item is trivial, since the sequence u{f) = ('Uj(/))j>i is bounded, say 
by the constant C. Indeed, then the Cesaro means of the sequence u(/)^ = 
(iij(/)^)j>i are bounded, for any ^ > by 

<T%Hff) < Ca%{Hf)\). (5.23) 
This implies the {C,9) convergence of the sequence u{f )'^. □ 

Thus, Lemma [5 . 2 1 shows that Theorem 12 . 31 implies Theorem II .11 Indeed, the 
general assumptions needed in Theorem 12.31 about the Cesaro convergence 
of the sequences u{f) and n(/)^ are satisfied by Lemma [5^21 

6. Proof of Theorem 12.41 and Theorem 11.51 

6.1. A simple Gaussian convergence result for series of Poisson 
random variables. 

We give here a result of convergence in distribution for the random variables 

N 

HN{f) = Y.W,u,{f) (6.1) 

i=i 

to a Gaussian law, once centered and normalized. This result is again el- 
ementary since it only uses the fact that the random variables Wj^s are 
independent and Poisson. 
Here, we assume as in Theorem 12.41 that 

oo 

Y,jRj{f? = ^ (6.2) 
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and that 

N 

i™<^ |uj| = ^max^ \jRj\ = oiC^jRjif ff^). (6.3) 

i=i 

Let us denote by 

N 

ri% = Var(/7jv) = 9Y,jRj{ff- (6-4) 

i=i 

Lemma 6.1. Under these assumptions, the distribution of 

hn - nHN] 

VN 

converges weakly to A/'(0, 1) as N ^ oo. 
Proof. Write i^^v for 

Hn - ^Hn] ^ Hn - E[gjv] ^ y u,{Wj - {9/j)) 

VN ^yVar[HN] jT[ Vn ' 

then 

N 



logE 



g.ti/ivj = ^_^(.itu,hN _ ituj/rjN - 1), (6.6) 



which gives the distribution of Hn in its Levy-Khintchine representation 
{aN,dMN,cr%), with 

«-|:(^(^)). (") 

MN = Y.-^u,/,r. (6-8) 

i=i 

and £7 AT = 0. 

We continue this proof by applying the Levy-Khintchine Convergence The- 
orem ([21]) P- 62): 

Consider a bounded continuous function / such that f{x) = for |x| < 6, 
then 

N 

H / 11 , \ 

(6.9) 

- .7 V 'IN I l^777K" 

J= 

Under the assumption ()6.3p . / fdMN = for large enough, so that 



[fdMN = Y.-f(^)l\u, 



hm / fdMN = I fdM = 0. (6.10) 

N~^ooJ J 



N-^oo _ 

Again, using the assumption ()6.3p . we have that for any ^ > 0, 



/ u- 
3 Vn 
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So for large enough, x^dM^ = 1 and so 



lim / x^dMN = 1. (6.12) 

N—^oo 



Moreover, for every define 



ej, = "^""^^^-^^ '"^ (6.13) 

VN 



then we can bound above by 



By the assumption (|6.3p . we see that 

hm a^v = 0. (6.15) 

By (I6TI)]1 . (I6J2I1 . (I6J5]1 and using Theorem 3.21, p. 62 in we see that 
H]\[ converges in distribution to the infinitely divisible distribution with 
Levy-Khintchine representation (a,M,a'^) = (0,0,1), i.e. to the standard 
normal Gaussian AA(0, 1). 

□ 

6.2. Proof of Theorem [23 

We proceed now to the proof of Theorem 12.41 by using the Feller coupling 
bounds proved in Section 4. Again, we assume here, as in Theorem l2.4l that 

oo 

Y,jR,iff = oo (6.16) 

i=i 

and that 

N 

max^|ji?,| = o((J] ji?,(/)2)5) (6.17) 



which can be rewritten 



max \uA = oiriM) (6.18) 
i<i<Af I J I \ ' / 



By Lemma 14.21 we know that 

E(\Gn - Hn\) <C(9) max |n,| = o(77Ar) (6.19) 

l<j<N ■' 

Again, denote by ■= and Gn := ^ilLi^ilil^ Then obviously, 

E(|G^-^;v)|) = o(l), (6.20) 
which, together with the convergence result (Lemma l6.ip for H]\f proves that 
Gn converges in distribution to a standard Gaussian law N(0,1). 

Moreover, this also proves that 

N 

e{Gn) = nHN)) + o{m) = eY,Rj + oim)- (6.21) 
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E{{Gn - Hn?) < C{e) ^max^'^i = o{r]l). (6.22) 

Then obviously, 

E((Gjv - Hn?) < 2(E((Gjv - Hn?) + E{Gn - Hn?) = o{rj%). (6.23) 
So that 



I v'Var(GAr) - \AM^I = oim) (6-24) 

and thus, 

Var(G7v) ~ Var(i/w) = ri%. (6.25) 
From these three results, we get that 

Gn - lE(Gjv) 
VVar(G^) 

converges in distribution to a standard Gaussian A/'(0, 1). But Ia,Nif) — 
IE(/(T,Ar(/)) has the same distribution as Gn — ^{Gn) and thus, 

/,,iv - IE(/<x,iv) 



converges also in distribution to a A/'(0, 1) distribution. We thus have proved 
the first statement of Theorem 12.41 Moreover, we have that 

E(/Ar) = / /(x)dx + E(GAr) = / /(x)dx + 6* V + o(r?Ar), (6.26) 
Jo Jo 

which is the second statement of Theorem 12.41 Finally, 

Var(/^,^) = Var(Gw) ~ Var(/7jv) = vl, (6-27) 

which is the third statement. We have completed the proof of Theorem 12.41 

6.3. Proof of Theorem II. 5L 

We prove here how Theorem 12.41 implies Theorem 11.51 In Theorem 11.51 we 
assumed that / is of bounded variation, which implies, as we have seen, that 
Uj{f) = 0(1), and thus, that 

max \uj\ = oiriM) (6.28) 
i<i<A^ \' / 

since the sequence t^at is assumed to diverge. This proves that the hypothesis 
of Theorem 11.51 are satisfied under those of Theorem 12.41 Thus we get that 
the conclusions of Theorem 12.41 are valid. They are almost exactly the same 
as the conclusions of Theorem 11.51 The only thing left to prove is the item 
(ii). But using Lemma [Ol 15.111 and the fact that Uj{f) = jRj{f) = 0{1), 
we have that 

„i N 

ni^Mf)]=N f{x)dx + y^Rjif) + 0{l). (6.29) 
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The bound 



N 



(6.30) 



Y,Rjif) = 0{logN) 

is trivial since again Rj = 0(i). With this we have derived Theorem 11.5 
from Theorem [2 



7. The expectation and the variance 

For the sake of completeness, we give here the explicit expressions for the 
expectation and the variance of Ia,N, when a is chosen from by the Ewens 
distribution with parameter 6. The basic computations for the expectation 
and the variance of the cycle counts can be simply derived by the following 
formula established by Watterson [35] (see Arratia, Barbour and Tavare [2], 
(4.7), p. 68): 

For every 6 > 1, (ri, . . . , r^) > 0, 



E 



N — m 



N 



-1 b 

n 

i=i 



(7.1) 



where m = Y^j=i ^^''^ ~ ^(^ — 1) . . . (x — r + 1) and 7 = — 1. 
Thus, the mean and the variance of I^^n can be easily computed. 

Lemma 7.1. 

N 



j-1 

nhAf)] = N fix)dx + y^na,{a)]uj{f) 

1 N 

= N / f{x)dx+ey2^N{j)Rjif) 
Jo ,=1 

(7.2) 

Proof. From the general formula ()7.ip . we easily see that for any j > 1 and 
any 6* > 0, 

ft 

(7.3) 



Thus, (|7.2p follows immediately. 

The variance of I^^n is given by the following lemma: 
Lemma 7.2. 



□ 



Var[/,,jv(/)] = Var 



N 



i=i 



i=i hj'<N 



R.Rf {"fNij + j')lj+f<N -'fN{j)'fN{j')) (7.4) 
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Proof. Again, from the general formula (jT.ip . we easily see that for any j > 1 
and any 6 > 0, 

02 

Mo^jO^j'] = JjJ^N^^ + f)h+j'<N (7.5) 

and 

Ee[a]] = -^^N{2j)lj<N/2 + -^Jv(j)lj<jv. (7.6) 
The variance of aj is therefore given by 

Yaieiaj] = -^N{j)lj<N + -2^N{2j)lj<^/2 - -2^N{jflj<N (7.7) 
and the covariance by 

02 q2 

Co-ve[aj,aj>\ = -j;^ n{3 + 3')lj+j'<N - —r,^ N{j)^ N{j')lj<Nlj'<N, (7.8) 

for j 7^ j' . Then, the variance of I^^n given in (j7.4p follows immediately. □ 
Remark 8. The case where 9 = 1 \s particularly simple. Indeed then 

1 

E[/,,jv(/)] = iV / f{x)dx + 9y^R,{f) 
Jo ~\ 
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and 

Var[/,,;v(/)] = V%- 

Thus the asymptotic formulae we give in this work are then exact. For 
general values of > 0, it is possible to derive these asymptotic expressions 
directly form the explicit formulae given in this Section, without using the 
bounds on the Feller coupling, but this is not a trivial matter, in particular 
when 9 < 1. 
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